Daily R - March 2021
Packages
20210301
rayshader
20210302
Tidy Tuesday - Super Bowl Commercials
20210303
more rayshader
20210304
Slacker
20210305
20210308
gghighlight
20210309
Tidy Tuesday - Bechdel Data
## Warning: Tried to calculate with group_by(), but the calculation failed.
## Falling back to ungrouped filter operation...## label_key: binary## Too many data points, skip labeling
## Warning in grid.Call(C_stringMetric, as.graphicsAnnot(x$label)): font family not
## found in Windows font database
## Warning in grid.Call(C_stringMetric, as.graphicsAnnot(x$label)): font family not
## found in Windows font database
## Warning in grid.Call(C_stringMetric, as.graphicsAnnot(x$label)): font family not
## found in Windows font database## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database
20210310
More Tidy Tuesday
## Warning: Outer names are only allowed for unnamed scalar atomic inputs
20210311
terra
## [1] "C:/Users/SimsCR/Documents/R/R-4.0.2/library/terra/ex/test.tif"## source nlyr
## 1 C:/Users/SimsCR/Documents/R/R-4.0.2/library/terra/ex/test.tif 1## [1] TRUE
20210312
20210315
India
20210316
Tidy Tuesday - Steam Games
20210317
Continue Tidy Tuesday
20210318
Open Street Maps Data. The plotting is slow.
20210319
terrainr and progressr
## Warning in runif(100, 47.5, 45.5): NAs produced20210322
Trying cverdel’s rayshader experiment
20210323
Tidy Tuesday - UN Votes
20210324
ggbeeswarm
20210325
ggforce and Palmer Penguins and hex colors
## `geom_smooth()` using formula 'y ~ x'
20210326
Modeling with Palmer Penguins using broom
##
## Call:
## lm(formula = bill_length_mm ~ species + sex + flipper_length_mm +
## island, data = .)
##
## Residuals:
## Min 1Q Median 3Q Max
## -6.6263 -1.2580 -0.0012 1.2035 12.1675
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 18.37374 4.07711 4.507 9.20e-06 ***
## speciesChinstrap 9.74463 0.42943 22.692 < 2e-16 ***
## speciesGentoo 5.67808 0.73776 7.696 1.68e-13 ***
## sexmale 3.00340 0.28938 10.379 < 2e-16 ***
## flipper_length_mm 0.10117 0.02191 4.618 5.59e-06 ***
## islandDream -0.59680 0.45688 -1.306 0.192
## islandTorgersen -0.18159 0.47721 -0.381 0.704
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 2.256 on 326 degrees of freedom
## (11 observations deleted due to missingness)
## Multiple R-squared: 0.8329, Adjusted R-squared: 0.8298
## F-statistic: 270.8 on 6 and 326 DF, p-value: < 2.2e-16## # A tibble: 7 x 5
## term estimate std.error statistic p.value
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 (Intercept) 18.4 4.08 4.51 9.20e- 6
## 2 speciesChinstrap 9.74 0.429 22.7 4.68e-69
## 3 speciesGentoo 5.68 0.738 7.70 1.68e-13
## 4 sexmale 3.00 0.289 10.4 5.41e-22
## 5 flipper_length_mm 0.101 0.0219 4.62 5.59e- 6
## 6 islandDream -0.597 0.457 -1.31 1.92e- 1
## 7 islandTorgersen -0.182 0.477 -0.381 7.04e- 1## GVIF Df GVIF^(1/(2*Df))
## species 12.849772 2 1.893319
## sex 1.369800 1 1.170384
## flipper_length_mm 6.151079 1 2.480137
## island 3.695930 2 1.386535## Analysis of Variance Table
##
## Response: bill_length_mm
## Df Sum Sq Mean Sq F value Pr(>F)
## species 2 7015.4 3507.7 689.2651 < 2.2e-16 ***
## sex 1 1135.7 1135.7 223.1630 < 2.2e-16 ***
## flipper_length_mm 1 109.4 109.4 21.4948 5.141e-06 ***
## island 2 9.4 4.7 0.9256 0.3973
## Residuals 326 1659.0 5.1
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Call:
## lm(formula = bill_length_mm ~ species + sex + body_mass_g + island,
## data = .)
##
## Residuals:
## Min 1Q Median 3Q Max
## -6.938 -1.474 0.002 1.258 10.490
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 31.2376758 1.3782006 22.666 < 2e-16 ***
## speciesChinstrap 10.2815103 0.4103907 25.053 < 2e-16 ***
## speciesGentoo 6.1490998 0.6724810 9.144 < 2e-16 ***
## sexmale 2.5327994 0.3613959 7.008 1.39e-11 ***
## body_mass_g 0.0017443 0.0003939 4.429 1.29e-05 ***
## islandDream -0.4635553 0.4574428 -1.013 0.312
## islandTorgersen 0.0922458 0.4744012 0.194 0.846
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 2.261 on 326 degrees of freedom
## (11 observations deleted due to missingness)
## Multiple R-squared: 0.8321, Adjusted R-squared: 0.829
## F-statistic: 269.2 on 6 and 326 DF, p-value: < 2.2e-16## # A tibble: 7 x 5
## term estimate std.error statistic p.value
## <chr> <dbl> <dbl> <dbl> <dbl>
## 1 (Intercept) 31.2 1.38 22.7 5.90e-69
## 2 speciesChinstrap 10.3 0.410 25.1 5.61e-78
## 3 speciesGentoo 6.15 0.672 9.14 6.65e-18
## 4 sexmale 2.53 0.361 7.01 1.39e-11
## 5 body_mass_g 0.00174 0.000394 4.43 1.29e- 5
## 6 islandDream -0.464 0.457 -1.01 3.12e- 1
## 7 islandTorgersen 0.0922 0.474 0.194 8.46e- 1## GVIF Df GVIF^(1/(2*Df))
## species 10.322775 2 1.792459
## sex 2.125880 1 1.458040
## body_mass_g 6.529361 1 2.555261
## island 3.634004 2 1.380691## Analysis of Variance Table
##
## Response: bill_length_mm
## Df Sum Sq Mean Sq F value Pr(>F)
## species 2 7015.4 3507.7 685.8748 < 2.2e-16 ***
## sex 1 1135.7 1135.7 222.0653 < 2.2e-16 ***
## body_mass_g 1 101.5 101.5 19.8389 1.159e-05 ***
## island 2 9.1 4.6 0.8944 0.4099
## Residuals 326 1667.2 5.1
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
20210329
penguins_raw and visdat
## # A tibble: 7 x 6
## effect group term estimate std.error statistic
## <chr> <chr> <chr> <dbl> <dbl> <dbl>
## 1 fixed <NA> (Intercept) 38.9 0.433 89.8
## 2 fixed <NA> speciesChinstrap penguin (Pyg~ 10.2 0.530 19.1
## 3 fixed <NA> speciesGentoo penguin (Pygosc~ 8.60 0.482 17.9
## 4 fixed <NA> islandDream -0.304 0.552 -0.552
## 5 fixed <NA> islandTorgersen 0.0409 0.586 0.0698
## 6 ran_pars individu~ sd__(Intercept) 1.50 NA NA
## 7 ran_pars Residual sd__Observation 2.55 NA NA## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: culmen_length_mm
## Chisq Df Pr(>Chisq)
## species 666.3898 2 <2e-16 ***
## island 0.4801 2 0.7866
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
20210330
Tidy Tuesday
## --- Compiling #TidyTuesday Information for 2021-03-30 ----## --- There are 5 files available ---## --- Starting Download ---##
## Downloading file 1 of 5: `ulta.csv`
## Downloading file 2 of 5: `sephora.csv`
## Downloading file 3 of 5: `allShades.csv`
## Downloading file 4 of 5: `allNumbers.csv`
## Downloading file 5 of 5: `allCategories.csv`## --- Download complete ---
20210331
Just working on completing the month/knitting
Ideas
terrainr, template, rvest, elevatr, gtsummary, officedown/officer, glue, beepr, guide_colorbar, tidymodels, ggsankey, gluedown,
babynames - https://towardsdatascience.com/how-popular-is-your-name-mini-data-viz-challenge-in-r-1dc31f6a22dc
https://dominicroye.github.io/en/2019/visualize-urban-growth/
https://github.com/Z3tt/OutlierConf2021/blob/main/R/OutlierConf2021_ggplotWizardry_HandsOn.Rmd
https://jokergoo.github.io/circlize_book/book/index.html
http://shadedrelief.com/SampleElevationModels/
https://www.neonscience.org/resources/learning-hub/tutorials
Extra
elevatr - Examples in Vignette don’t work.